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Abstract 

We show that a large class of dissipative systems can be brought to a 
canonical form by introducing complex co-ordinates in phase space and a 
complex-valued hamiltonian. A naive canonical quantization of these sys- 
tems lead to non-hermitean hamiltonian operators. The excited states are 
unstable and decay to the ground state . We also compute the tunneling 
amplitude across a potential barrier. 
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1 Introduction 



In many physical situations, loss of energy of the system under study to the out- 
side environment cannot be ignored. Often, the long time behavior of the system 
is determined by this loss of energy, leading to interesting phenomena such as 
attractors. 

There is an extensive literature on dissipative systems at both the classical and 
quantum levels (See for example the textbooks flU [2l (3J ) . Often the theory is based 
on an evolution equation of the density matrix of a 'small system' coupled to a 
'reservoir' with a large number of degrees of freedom, after the reservoir has been 
averaged out. In such approaches the system is described by a mixed state rather 
than a pure state: in quantum mechanics by a density instead of a wavefunction 
and in classical mechanics by a density function rather than a point in the phase 
space. 

There are other approaches that do deal with the evolution equations of a pure 
state. The canonical formulation of classical mechanics does not apply in a direct 
way to dissipative systems because the hamiltonian usually has the meaning of 
energy and would be conserved. By redefining the Poisson brackets flU , or by 
using time dependent hamiltonians J51, it is possible to bring such systems within 
a canonical framework. Also, there are generalizations of the Poisson bracket that 
may not be anti-symmetric and/or may not satisfy the Jacobi identity (6l|71 which 
give dissipative equations. 

We will follow another route, which turns out in many cases to be simpler than 
the above. It is suggested by the simplest example, that of the damped simple 
harmonic oscillator. As is well known, the effect of damping is to replace the 
natural frequency of oscillation by a complex number, the imaginary part of which 
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determines the rate of exponential decay of energy. Any initial state will decay to 
the ground state (of zero energy) as time tends to infinity. The corresponding co- 
ordinates in phase space (normal modes) are complex as well. This suggests that 
the equations are of hamiltonian form, but with a complex- valued hamiltonian. 

It is not difficult to verify that this is true directly. The real part of the hamil- 
tonian is a harmonic oscillator, although with a shifted frequency; the imaginary 
part is its constant multiple. If we pass to the quantum theory in the usual way, 
we get a non-hermitean hamiltonian operator. Its eigenvalues are complex valued, 
except for the ground state which can be chosen to have a real eigenvalue. Thus 
all states except the ground state are unstable. Any state decays to its projection to 
the ground state as time tends to infinity. This is a reasonable quantum analogue 
of the classical decay of energy. 

We will show that a wide class of dissipative systems can be brought to such 
a canonical form using a complex-valued hamiltonian. The usual equations of 
motion determined by a hamiltonian and Poisson bracket are 

d ( P\ ( i H >P} 



dt 



x J \ {H, x} 



(1) 



At first a complex- valued function 

H = H 1 + iH 2 (2) 
does not seem to make sense when put into the above formula: 

dt\ x ) [{H^x}) \{H 2 ,x}) 

since the l.h.s. has real components. How can we make sense of multiplication by 
% and still get a vector with only real components? 







f-v\ 









Let us consider a complex number z = x + iy as an ordered pair of real 
numbers (x, y). The effect of multiplying z by z is the linear transformation 



(4) 



on its components. That is, multiplication by i is equivalent to the action by the 
matrix 

Note that J 2 = —1. Geometrically, this corresponds to a rotation by ninety de- 
grees. 

Generalizing this, we can interpret multiplication by i of a vector field in phase 
space to mean the action by some matrix J satisfying 

J 2 = -1. (6) 

Given such a matrix, we can define the equations of motion generated by a complex- 
valued function H = Hi + iH 2 to be 

dt{ x ) \{H u x}) + \{H 2 ,x}) 

Our point is that the infinitesimal time evolution of a wide class of mechanical 
systems is of this type for an appropriate choice of {, }, J, Hi and H 2 . 

In most cases there is a complex-cordinate system in which J reduces to a 
simple multiplication by i; for example on the plane this is just z = x + iy. For 
such a co-ordinate system to exist the tensor field has to satisfy certain integrability 
conditions in addition to © above. These conditions are automatically satisfied if 
the matrix elements of J are constants. 
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What would be the advantage of fitting dissipative systems into such a com- 
plex canonical formalism? A practical advantage is that they can lead to better 
numerical approximations, generalizing the symplectic integrators widely used 
in hamiltonian systems: these integrators preserve the geometric structure of the 
underlying physical system. Another is that it allows us to use ideas from hamil- 
tonian mechanics to study structures unique to dissipative systems such as strange 
attractors. We will not pursue these ideas in this paper. 

Instead we will look into the canonical quantization of dissipative systems. 
The usual correspondence principle leads to a non-hermitean hamiltonian. As in 
the elementary example of the damped simple harmonic oscillator, the eigenvalues 
are complex- valued. The excited states are unstable and decay to the ground state. 
Non-hermitean hamiltonians have arised already in several dissipative systems in 
condensed matter physics [9] and in particle physics ifTOl . The Wigner-Weisskopf 
approximation provides a physical justification for using a non-hermitean hamil- 
tonian. A dissipative system is modelled by coupling it to some other 'external' 
degrees of freedom so that the total hamiltonian is hermitean and is conserved. In 
second order perturbation theory we can eliminate the external degrees of freedom 
to get an effective hamiltonian that is non-hermitean. 

It is interesting to compare our approach with the tradition of Caldeira and 
LeggettjU. Dissipation is modelled by coupling the original ('small') system to 
a thermal bath of harmonic oscillators. After integrating out the oscillators in the 
path integral formalism an effective action for the small system is obtained. A 
complication is that this effective action is non-local: its extremum ( which dom- 
inates tunneling) is the solution of an integro-differential equation. We will see 
that the integral operator appearing here is also a complex structure (the Hilbert 
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transform), although one non-local in time and hence different from our use of 
complex structures. 

We calculate the tunneling amplitude of a simple one dimensional quantum 
system within our framework. Dissipation can increase the tunneling probability, 
which is not allowed in the Caldeira-Leggett model. 

We begin with a brief review of the most elementary case, the damped simple 
harmonic oscillator. Then we generalize to the case of a generic one dimensional 
system with a dissipative force proportional to velocity. Further generalization to 
systems with several degrees of freedom is shown to be possible provided that the 
dissipative force is of the form 

- d a d b W^- (8) 
at 

for some function W. In simple cases this function is just the square of the dis- 
tance from the stable equilibrium point. Finally, we show how to bring a dissipa- 
tive system whose configuration space is a Riemannian manifold into this frame- 
work. This is important to include interesting systems such as the rigid body or 
a particle moving on a curved surface. We hope to return to these examples in a 
later paper. 

2 Dissipative Simple Harmonic Oscillator 

We start by recalling the most elementary example of a classical dissipative sys- 
tem, described the differential equation 

x + 27a; + u 2 x = 0, 7 > 0. (9) 
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We will consider the under-damped case 7 < lu so that the system is still oscilla- 
tory. 

We can write these equations in phase space 



x = v (10) 
p = —2^p — u 2 x (11) 



H = \[p 2 + u 2 x 2 } (12) 



The energy 

u — _ 

2 

decreases monotonically along the trajectory: 

dH 

— = pp + uj 2 xx = -2>yp 2 < 0. (13) 
dt 

The only trajectory which conserves energy is the one with p — 0, which must 
have x = as well to satisfy the equations of motion. 

These equations can be brought to diagonal form by a linear transformation: 

dz 

z = A \-i(p + 7X) + uix] , — = [—i + uoi]z (14) 

where 



u; 1 = v ^ 2 -7 2 . (15) 

The constant A that can be chosen later for convenience. These complex co- 
ordinates are the natural variables (normal modes) of the system. 

2.1 Complex Hamiltonian 

We can think of the DSHO as a generalized hamiltonian system with a complex- 
valued hamiltonian. 
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The Poisson bracket {p, x} — 1 becomes, in terms of the variable z, 

{z*,z} = 2iuj 1 \A\ 2 (16) 

So if we choose A = - 7 = 

{z*,z} = i (17) 

So the complex-valued function 

H = (u)i + i-y)zz*. (18) 

satisfies 

§ = <"•*>• f = <"••*•> (19) 

Of course, the limit 7^0 this tends to the usual hamiltonian H = uzz*. 
Thus, on any analytic function ip, we will have 

^ = {H^}=[u Jl + t 1 ]z^ (20) 

2.2 Quantization 

By the usual rules of canonical quantization, the quantum theory is given by turn- 
ing H into a non-hermitean operator by replacing z 1— > at, z* ^ Tta and 

[a, a*] = 1, a* = z, a = — , Ti = %{uji + i^)a} a. (21) 

The effective hamiltonian H = Hi + iH 2 is normal ( i.e., its hermitean and anti- 
hermitean parts commute, [H 1 ,H 2 ] = ) so it is still meaningful to speak of 
eigenvectors of H. The eigenvalues are complex 

(uj l + i 1 )n,n = 0,1,2,- ••. (22) 
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The higher excited states are more and more unstable. But the ground state is 
stable, as its eigenvalue is zero. 
Thus a generic state 

oo 

= (23) 

n=0 

will evolve in time as 

oo 

^(t) = ]T ifj n e lhluJ+ ^ ]nt \n > . (24) 

n=0 

Unless ip happens to be orthogonal to the ground state |0 >, the wavefunction will 
tend to the ground state as time tends to infinity; final state will be the projection 
of the initial state to the ground state. This is the quantum analogue of the classical 
fact that the system will decay to the minimum energy state as time goes to infinity. 
All this sounds physically reasonable. 



2.3 The Schrodinger Representation 



In the Schrodinger representation,this amounts to 



1 



a = 



y/2hu>i 
H = 



* 9 

ujix + n— 

ox 



I 9 

uj\x — n— 
ox 



l + i 



jL 



h 2 d 2 

12 2 1 + 



Thus the operator representing momentum p is 

p = —in— 7X 

Ox 

which includes a subtle correction dependent on the friction. 
The time evolution operator can be chosen to be 



(25) 
(26) 

(27) 



T^Schr = H + H, 



diss 



(28) 



where 

h 2 d 2 



H =-T^ + 1 ^ (29) 



is the usual harmonic oscillator hamiltonian and 

h 2 d 2 



ZJ 1 2 2 i • 7 



2 dx 2 



+ \uo\x 2 - \TlUJl 



(30) 



This is slightly different from the operator H above, because the ground state 
energy is not fixed to be zero. The constant in H diss has been chosen so that this 
state has zero imaginary part for its eigenvalue. 

3 Dissipative System of One Degree of Freedom 

We will now generalize to a non-linear one-dimensional oscillator with friction: 

dp dV dx „ ,„, N 

The DSHO is the special case V(x) = \u?x 2 . The idea is that we lose energy 
whenever the system is moving, at a rate proportional to its velocity. It again 
follows that 

IT = - 2lp2 ~ (32) 



where H = \p 2 + V. These equations can be written as 



-± = {H,?}-i»d j H (33) 



2 7 0' 

where 7=1 I is a positive but degenerate matrix. 
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3.1 Complex Effective Hamiltonian 

In the case of the DSHO, we saw that it is the combination 

Pi = P + ix. (34) 

rather than p that appears naturally ( for example in the normal co-ordinate z). In 

terms of (pi , x) the equations above take the form 

dpi dVi dx 

_ 7Pi> -77= Pi ~l x i (35) 



dt dx " dt 



where 



i.e. 



V x {x) = V(x) - Wx 2 - (36) 



We would like to see if we can write these as canonical equations of motion 
with a complex- valued effective hamiltonian. Suppose we define 



tf^-L^V], 7=1 "I (38) 

for some constant uj 2 which we will choose later. Note that J is a complex struc- 
ture; i.e., 

j2 = ~(l J- (39) 

Then the equations of motion take the form 

d (Pi\ ( {Hi,pi}\ ( {H 2 ,pi}\ 

J I • (40) 



-u 2 



dt 



x ) \ {Hi,x} J \ {H 2 ,x} 



In terms of the complex co-ordinate 
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this becomes 

^ = {H 1+i H 2 ,z}. (42) 

Thus the non-linear oscillator also can be written as a canonical system with a 
complex valued hamiltonian H = Hi + %H 2 \ with {Hi, H 2 } ^ in general. But 
there are many ways to do this, parametrized by uo 2 . The natural choice is 

U2 = yV"(xo) — 7 2 > (43) 

where x is the equilibrium point at which V'(x ) = 0. Then, in the neighborhood 
of the equilibrium point the complex structure reduces to that of the DSHO. 

3.2 Quantization of a Dissipative One-dimensional System 

We can quantize the above system by applying the usual rules of canonical quan- 
tization to the complex effective hamiltonian Hi + iH 2 to get the operator: 



fr 2 ff 2 -v 
H = -~ + V(x)-tfa? + il- 
2 ox z uj 2 

The Schrodinger equation 



k 2 d 2 , 2 2 

+ 7jUJ 2 X + c 



2 dx 2 



(44) 



- ih^- = Hi) (45) 
at 

then determines the time evolution of the dissipative system. 

The anti-hermitean part of the hamiltonian is bounded below so the imaginary 
part of the eigenvalues will be bounded below. We can choose the real constant 
c such that the eigenvalue with the smallest imaginary part is actually real. Then 
the generic state will evolve to this ' ground state'. 
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Remarks 



1 . Our model of dissipation amounts to adding a term 

n 1 2 2^-7 r h 2 d 2 



^2 

to the usual conservative hamiltonian 



2 dx 2 



12 2 



+ c (46) 



* = -f^ + ^- (47) 

The dissipative term we add is very close to being anti-hermitean; i.e., ex- 
cept for the term — \^ 2 x 2 which is second order in the dissipation. 

2. The Schrodinger equation for time reversed QM ( evolving to the past rather 
than future) is H) . The eigenvalues of 1~0 are the complex conjugates of 
those of TC, but the eigenfunctions may not be the same as those of H since 
[H,Tft] 7^ in general; i.e., the operator H may not be normal. ( For 
the DSHO the commutator vanished so the issue did not arise.) Moreover 
the eigenfunctions of H corresponding to distinct eigenvalues need not be 
orthogonal; they would still be linearly independent, of course. 

There are examples of such non-normal hamiltonians in nature when time- 
reversal invariance is violated ( e.g., KK oscillation)[10|. But in ordinary 
quantum mechanics such loss of time reversal invariance might be unset- 
tling. Carl Bender [fTTTl has suggested that for such non-hermitean hamil- 
tonians, real eigenvalues and time reversal invariant dynamics can be re- 
covered by modifying the inner product in the Hilbert space. This is the 
quantum counterpart to the reformulation of classical DSHO as a conser- 
vative canonical system by modifying the Poisson bracket and hamiltonian 
0. 
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However we note that in the presence of dissipation, the classical equations 
of motion are not time reversal invariant either: energy would grow rather 
than dissipate. So we should not expect quantum mechanics of dissipative 
systems to be time reversal invariant either. The appropriate symmetry is 

[ft( T )]t = n{-y) (48) 

which is satisfied in our case. 

3. Note that within our model, the anti-hermitean part of the hamiltonian is a 
sort of harmonic oscillator even if the hermitean part has non-linear classical 
dynamics. This is because we chose a particularly simple form of dissipa- 
tion, p ~ — 7p. If we had chosen a more complicated ( e.g., non-linear) 
form of dissipation, the anti-hermitean part would be more complicated. It 
is common to model a dissipative system by coupling it to a thermal bath of 
oscillators. The dissipation is determined by the spectral density of the fre- 
quencies of these oscillators. Each choice of spectral density leads to a dif- 
ferent dissipation term and, in our description, to a different anti-hermitean 
part for the hamiltonian. But if the dissipative force is small it is reasonable 
to expect that it is linear in the velocity. 

4 Comparison with the Approach of Caldeira and 
Leggett 

It is interesting to see how our approach compares with the influential work of 
Caldeira and Leggett [8] 
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Their basic idea is to couple the one dimensional system to a collection of har- 
monic oscillators; the transfer of energy to the oscillators leads to its dissipation 
in the original system. The lagrangian is3 

L=\x 2 - V(x) + Y\\m a ql - > a c^] - x£C Q g Q - \/\u 2 x 2 . (49) 

a a 

In the last term, Ao; 2 is chosen such that the equilibrium point of the variable x is 
still the minimum of V: it is a finite renormalization of the potential induced by 
the coupling to the oscillators. 

In the euclidean path integral approach, they derive an effective action after 
integrating out these oscillators: 

rT /*oo rT 

S eS [x}= [\x 2 + V{x)}dr + ± / dr' dT[x(r)-x{T')} 2 a(r-T') (50) 

JO J-oo JO 



where 

CI 
a 4m a u a 

All the properties of the oscillators that matter to us are contained in the function 



«W = E -^e— = - / p(u)e-^du (51) 



p(u). The choice p(u) = corresponds to the usual dissipative term . This 
at first sounds strange: how can this effective action that is non-local in time be 
related to the differential equation of the damped oscillator? The minimum of S e g 
( which determines the tunneling amplitude in the leading WKB approximation) 
is given by the integro-differential equation^ 

<Px = 8V | 7p d7 j x(r) - x(t') 
dr 2 dx 7r 7-oo (r — r') 2 

In the absence of damping we would get this equation for the most likely path in 
the classically forbidden region by the replacement |fT2l : 

t ^ it (53) 



1 We change notation slightly from them to avoid conflict with ours. 

2 Here, V stands for the Cauchy Principal part. 
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in Newton's law: 

^ = J™. (54) 

dt 2 dx 

Such a naive replacement would not make sense in the damped equation: 

cPx dV dx d 2 x dV . dx ^ 
dt 2 dx dt dr 2 dx dr 

Since x is real, the multiplication by i on the r.h.s. does not make sense. 

Perhaps multiplication by i should be replaced by the action of some linear 
operator on the space of functions of one variable? If so the square of this lin- 
ear operator must be —1. There is a well-known operator with this property, the 
Hilbert transform [fT3l : 

Jx (t) = -V r ^ldt> (56) 

7T J-oo t — t 

The Caldeira-Legget equation of motion is obtained by replacing the % in the naive 
equation above by the operator J: 

d 2 x dV d g 
dr 2 dx ^ dr 

Thus although obtained by very different arguments, there is a way to under- 
stand the Caldeira-Legget equation of motion in terms of a complex structure ( the 
Hilbert transform) on the space of paths. We instead introduce a complex structure 
in the phase space at the classical level. Could ours be the canonical quantization 
of their non-local path integral? 



5 Tunneling in a Dissipative System 

Perhaps the most interesting question about quantum dissipative systems is how 
dissipation affects tunneling. The standard WKB approximation method adapts 
easily to our case. 

16 



For illustrative purposes, it suffices to consider a one-dimensional system with 
a potential 

V(x) = \u?x 2 (l - ^ (58) 

We ask for the tunneling probability amplitude from the origin x = to the point 
x = a in a long time. In the absence of dissipation this is given in the WKB 
approximation by 

e -i/;^. (59) 
The integral in the exponent is the minimum of the imaginary time action 

J [\x 2 + V(x(r))]dr (60) 

among all paths satisfying the boundary conditions x(0) = and x(oo) = a. 
(This minimizing path is the 'instanton'.) 

If we apply the WKB approximation to the Schrdinger equation we get 



Solving this, 



dx 



002 



1 1 9( ^\ J_ 1 2 

Ate) +w 



The tunneling probability is given by 



2Vi(x) + ryw 2 £ 2 



dx 



0. (61) 



(62) 



3 -2Re <f>(b) 



(63) 



where b is the point of escape from the potential barrier; it might depend on the 
dissipation. 



It looks most natural to choose o> 2 = = V u2 ~ 7 2 as m the case of the 
DSHO. Then, 



(j>{x) 



W X 
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xdx. 



a 



(64) 



Now, 



/ x\/b — xdx = — b 2 
o 15 



(65) 



There are small discrepancies ( up to higher order terms in 7) depending on 

LO 2 

whether we think of the point of escape as a, or -^a or the complex number 



This differs from the results of Caldeira and Legget: the tunneling probability is 
enhanced by dissipation. There could still be systems in nature that are described 
by our model. 

6 Multi-Dimensional Dissipative Systems 

Now we will generalize to a multi-dimensional dynamical system while also al- 
lowing for a certain kind of non-linearity in the frictional force. 

Suppose the hamiltonian is the sum of a kinetic and a potential energy, 




H= \p a Pa + V{x). 



(67) 



With the usual Poisson brackets 



{ Pa ,x b } = 5 b a , { PaiPb } = = { x ^ x b } 



(68) 



the conservative equations of motion are 




dx' 



dt 



= P 



(69) 
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We now assume that the frictional force is of the form —2j ab x b for some pos- 
itive non-degenerate tensor ^ ab that might depend on x. The idea is that the sys- 
tem loses energy when parts of it move relative to other parts or relative to some 
medium in which the system is immersed. Then the equations of motion become 

^ = -d a V-2 labPb , d -^=V a - (70) 

Here we raise and lower indices using the flat euclidean metric S ab . So 

^ = -2 7 a W < 0. (71) 

If the dissipation tensor happens to be the Hessian of some convex function 

lab = d a d b W, (72) 

it is possible to write these as Hamilton's equations with a complex-valued hamil- 
tonian. 

Then, we would have 

j t d a W = lab p\ lab d b W = d a (\d b Wd b W). (73) 

Using these identities we can rewrite the equations of motion in the new variables: 

Pa=Pa + a W (74) 



as 

dx a 

dt ' J ~dt 



- = -d a [V - \{dWf] - lab p b) ^= Pa -d a W. (75) 



This motivates us to define 



Hl = ±f + V- \(dW)\ H 2 = —[\ lab p a p b + u 2 2 W], J 



(76) 
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for some positive constant u 2 . Note that J is a complex structure, J 2 = — 1. Then 

d(p\ f{Hi,P}\ T ({H2,P} i 

+ Jl I (77) 



dt 



xj \{H 2 ,x} 



In terms of the complex variable 
this is just 

^ = {H 1 + iH 2 ,z a } (79) 

just as before. Thus once the dissipation tensor is of the form = d a dbW the 
whole framework generalizes easily. The effect of dissipation is to add the term 

tfdiss = -\{dW? + —Ih-bPaPb + u 2 2 W] (80) 

L0 2 

to the hamiltonian. The classical equations turn out to be independent of the 
choice of uj 2 , but the quantum theory will depend on its choice. 



6.1 Quantization 

We can then quantize this operator in the Schrodinger picture: 

7i = H + //diss, (81) 

where 



U! 2 



(82) 

A technical point to note here is that there are two competing metrics in the 
story. The kinetic energy is \p a p a , determined by the euclidean metric. But the 
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quadratic term in the dissipative part of the hamiltonian is determined by some 
other tensor ^ ab . We have chosen to use the Euclidean metric 5 ab to define deriva- 
tives and to raise and lower indices. Thus, r y ab = 5 ac 5 bd ^ c d and not the inverse of 

lab- 

The operator d a (i ab d b if?j is thus a kind of 'mixed' laplacian that uses j ab as 
well as, implictly, the euclidean metric. Since d a ^ bc 7^ in general, the ordering 
of factors is important, with the order we chose,<9 a (l ab d b ip) is hermitean and 
positive. 

7 Geometric Model of Dissipative Mechanics 

We will now further generalize to the case of dissipative dynamics on a cotangent 
bundle T*Q. 

The hamiltonian is again the sum of kinetic and potential energies 

H = \g ab p aPb + V(x). (83) 

except that we now allow the tensor g ab not to be constant. We will use g ab ( the 
inverse of g ab ) as the Riemann metric, used to define covariant derivative V a and 
to raise and lower indices. 

The conservative equations of motion are 

^ + r a bc p b p c = -g ab d b V, ^=P a (84) 

where Tf c is the usual Christoffel symbol. 
With friction added, 

^ + r a bcPb p c = -g ab d b V - 2 7 a V ^ = P a (85) 
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where ^ ab is the dissipation tensor, assumed to be positive and non-degenerate. 
Again, if the dissipation tensor is the Hessian of a convex function 

lab = 9ac9bdl bd = V a d b W, (86) 
there are simplifications because 

j t (V a W) + T a bc V 6 Wp c = labV b (87) 



We define again 



to get 



where 





(88) 



(89) 



Hi = \g ab VaVb + V-{VW)\ H 2 = -L^PaPb + ^W, J = 

ZiOo 




(90) 

Since again J 2 = — it is an almost complex structure; but it may not be integrable 
in general. Every tangent bundle has a natural almost complex structure; J is 
simply its translation to the co-tangent bundle T*Q using the metric g ab which 
identifies the tangent and co-tangent bundles. 

Because J may not be integrable, we are not able to rewrite this in terms 
of a complex co-ordinate z in general. Nevertheless, we can think of the above 
equations as a generalization of hamilton's equations to a complex hamiltonian 
#i + iH 2 . 

Clearly, it is possible to quantize these system by applying the correspondence 
principle. Since the ideas are not very different, we will not work out the details. 
The hamiltonian is: 
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H = H + H diss , (91) 

where 



U>2 



(92) 

The imaginary constant ic is chosen such that the ground state is stable. 
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